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The  results  from  a  —  1  km  resolution  HYbrid  Coordinate  Ocean  Model  (HYCOM).  forced  by  1/2  Navy 
Operational  Global  Atmospheric  Prediction  System  (NOGAPS)  atmospheric  data,  were  used  in  order  to 
study  the  dynamic  response  of  the  Persian  Gulf  to  wintertime  shamal  forcing.  Shamal  winds  are  strong 
northwesterly  winds  that  occur  in  the  Persian  Gulf  area  behind  southeast  moving  cold  fronts.  The 
period  from  20  November  to  5  December  2004  included  a  well  defined  shamal  event  that  lasted 
4-5  days.  In  addition  to  strong  winds  (16  m  s_1)  the  winter  shamal  also  brought  cold  dry  air  (T„= 20  C. 
qa=  10  g  kg  ’)  which  led  to  a  net  heat  loss  in  excess  of  1000  W  m“ 2  by  increasing  the  latent  heat  flux. 
This  resulted  in  SST  cooling  of  up  to  10  X  most  notably  in  the  northern  and  shallower  shelf  regions. 
A  sensitivity  experiment  with  a  constant  specific  humidity  of  qa=  15  g  kg  1  confirmed  that  about  38%  of 
net  heat  loss  was  due  to  the  air-sea  humidity  differences.  The  time  integral  of  SST  cooling  closely 
followed  the  air-sea  heat  loss,  indicating  an  approximate  one-dimensional  vertical  heat  balance.  It  was 
found  that  the  shamal  induced  convective  vertical  mixing  provided  a  direct  mechanism  for  the  erosion 
of  stratification  and  deepening  of  the  mixed  layer  by  30  m.  The  strong  wind  not  only  strengthened  the 
circulation  in  the  entire  Persian  Gulf  but  also  established  a  northwestward  flowing  Iranian  Coastal 
Current  (ICC.  25-30  cm  s”1)  from  the  Strait  of  Hormuz  to  about  52°E,  where  it  veered  offshore.  The 
strongest  negative  sea  level  of  25-40  cm  was  generated  in  the  northernmost  portion  of  the  Gulf  while 
the  wind  setup  against  the  coast  of  the  United  Arab  Emirates  established  a  positive  sea  level  of 
15-30  cm.  The  transport  through  the  Strait  of  Hormuz  at  56.2  E  indicated  an  enhanced  outflow  of  0.25 
Sv  (Sv*  106  m1  s  ')  during  24  November  followed  by  an  equivalent  inflow  on  the  next  day. 

c  2010  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

The  Persian  Gulf  (also  known  as  the  Arabian  Gulf)  is  a  shallow, 
semi-enclosed  basin  with  a  mean  depth  of  about  35  m  and  is 
connected  to  the  deep  (>  1000  m)  Gulf  of  Oman  through  the 
narrow  Strait  of  Hormuz.  It  consists  of  a  shallow  and  wide  shelf 
along  the  southern  part  of  the  Gulf  with  a  depth  less  than  50  m, 
and  a  deep  northern  part  with  a  maximum  depth  of  100  m 
separated  from  the  coast  by  a  narrow  shelf.  The  circulation  in  the 
Persian  Gulf  is  driven  by  five  types  of  forcing:  wind-stress,  surface 
buoyancy  fluxes,  freshwater  runoff,  water  exchange  through  the 
Strait  of  Hormuz  and  tides.  The  winds  in  the  Persian  Gulf  are 
predominantly  northwesterly  throughout  the  year.  During  winter 
(November-February)  the  winds  are  slightly  stronger (~5  ms1) 
than  those  during  the  summer  (June- September)  (~3  m  s  ’J.The 
cold,  dry-air  emanating  from  the  northwest  during  winter  leads  to 
excessive  evaporation  which  results  in  the  convective  formation 
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of  the  saltier  and  the  denser  Persian  Gulf  Water  (PGW).  The  tidal 
influence  on  the  circulation  in  the  Persian  Gulf  is  insignificant 
except  in  the  Strait  of  Hormuz  and  along  the  Iranian  coast  (Blain, 
1998,  2000). 

The  general  circulation  of  the  Persian  Gulf  is  cyclonic,  which  is 
bounded  by  an  Iranian  Coastal  Current  (ICC)  along  the  northern 
side,  flowing  northwestward  from  the  Strait  of  Hormuz  with 
speed  greater  than  lOcms-1  (Hunter,  1983),  and  a  south¬ 
eastward  current  in  the  southern  part  of  the  Gulf  (Figure  34, 
Reynolds.  1993).  The  ICC,  which  flows  against  the  prevailing 
northwesterly  winds,  is  primarily  driven  by  the  pressure  gradient 
(Chao  et  al.,  1992;  Lardner  et  al.,  1993;  Blain.  2000;  Thoppil  and 
Hogan,  submitted  for  publication).  The  most  important  circulation 
feature  is  the  transformation  of  ICC  into  a  series  of  mostly  cyclonic 
eddies  during  August-September  (the  Iranian  Coastal  Eddies,  ICE, 
Thoppil  and  Hogan,  submitted  for  publication).  Baroclinic  in¬ 
stability  along  the  ICC  appears  to  be  the  mechanism  behind  the 
production  of  these  eddies. 

Circulation  and  water  mass  formation  in  the  Persian  Gulf  have 
been  the  subject  of  several  previous  modeling  studies  (Lardner 
et  al.,  1993;  Chao  et  al.,  1992;  Blain,  2000;  Sadrinasab  and  K3mpf, 
2004;  K3mpf  and  Sadrinasab,  2006;  Azam  et  al.,  2006).  Several  of 
these  studies  used  hydrodynamic  numerical  models  or  spectral 
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models  of  coarse  to  moderate  horizontal  resolutions  (20-5  km), 
and  were  driven  by  monthly  mean  density  or  atmospheric 
forcing.  While  these  models  reproduced  many  characteristics  of 
circulation  and  water-mass  reasonably  well,  they  did  not  resolve 
the  short-term  changes  associated  with  the  synoptic  weather 
events.  Although  the  prevailing  winds  are  northwesterly 
throughout  the  year,  they  are  occasionally  interrupted  by  the 
episodic,  short-lived,  strong  synoptic  wind  events,  called  the 
shamal.  These  events  can  cause  abrupt  changes  in  the  circulation 
and  heat-budget  for  a  short  period  of  time  as  demonstrated  in  the 
Adriatic  Sea  responding  to  the  Bora  wind  events  (Lee  et  al.,  2005). 
None  of  the  above  studies  have  examined  the  response  of  the 
short-term  shamal  winds  on  the  circulation  and  heat-budget  in 
the  Gulf. 

The  goal  is  to  study  the  dynamic  responses  of  the  Persian  Gulf 
to  a  sudden  outbreak  of  a  shamal  event  with  particular  emphasis 
on  circulation  and  heat-budget.  The  approach  is  to  use  results 
from  a  ~1  km  horizontal  resolution  Hybrid  Coordinate  Ocean 
Model  (HYCOM)  of  the  Persian  Gulf.  The  rest  of  the  paper  is 
organized  as  follows:  Shamal  winds  are  introduced  in  Section  2. 
Section  3  describes  the  model  configuration  and  forcing  fields. 
Responses  of  surface  heat  fluxes,  sea  surface  temperatures  (SSTs), 
circulation  and  sea  surface  height  associated  with  the  shamal 
event  are  discussed  in  Section  4.  Finally,  Section  5  offers  a 
conclusion  of  the  study. 


2.  The  shamal  wind  event 

The  shamal  (Arabic  word  for  north)  wind  is  a  strong  northerly 
or  northwesterly  wind  event  that  occurs  in  the  Persian  Gulf 
during  winter  as  well  as  summer  (Perrone,  1979).  The  winter 
shamal  occurs  mainly  between  November  and  March  and  is 
directly  associated  with  mid-latitude  disturbances  moving  from 
west  to  east.  They  are  caused  by  the  presence  of  a  large  pressure 
gradient  that  develops  behind  a  cold  front  passage.  Upper  level 
subsidence  and  rapidly  building  high  pressure  over  Saudi  Arabia 
and  Iraq  reinforces  the  low-level  northwesterly  winds.  Based  on 
duration,  there  are  two  types  of  winter  shamal:  those  which  last 
24-36  h  and  those  which  last  for  a  longer  period  of  3-5  days. 
During  a  shamal,  winds  typically  reach  15-20  ms  1  near  the 
surface  during  the  daytime.  A  shamal  occurs  first  in  the 
northwestern  part  of  the  Gulf  and  then  spreads  southeast  behind 
the  advancing  cold  front.  Moderate  to  strong  winds  can  raise 
desert  surface  material  and  reduce  visibility.  Shamal  wind 
generated  storm  surges,  coupled  with  tidal  effects  can  lead  to 
significant  changes  in  the  sea  level  of  several  meters  (El-Sabh  and 
Murty.  1989). 


In  this  study,  the  period  from  20  November  to  5  December 
2004  was  selected  because  it  included  a  well  defined  shamal 
event  that  lasted  4-5  days  (24-29  November).  Fig.  1  shows  the 
cross-calibrated,  multi-platform  (CCMP,  Atlas  et  al.,  2009)  surface 
winds  illustrating  the  evolution  of  shamal  event.  The  CCMP 
wind  is  produced  by  combining  all  satellite  derived  wind-speed 
and  vector  observations,  combined  with  ship  and  buoy  data  and 
the  European  Centre  for  Medium-Range  Weather  Forecasts 
(ECMWF)  analyses  (see  Atlas  et  al.,  2009  for  details).  On  20 
November  2004,  prior  to  the  shamal  event,  the  winds 
were  northwesterly  or  southeasterly  with  wind-speeds  under 
Sms  \  When  the  shamal  began  to  blow  on  23  November,  the 
wind  direction  changed  from  southeasterly  to  a  strong 
northwesterly  and  on  24  November  12Z  the  wind-speed 
exceeded  15  ms  1.  The  maximum  wind-speed  (10-15  ms  }) 
occurred  in  the  alongshore  direction  away  from  the  coast  due  to 
orographic  control  on  either  side  of  the  Gulf.  The  wind  effect  was 
confined  to  the  region  west  of  54  E  (Fig.  1)  with  a  maximum 
occurring  north  of  Qatar. 

A  comparison  of  time-series  of  wind-stress  (N  m  2)  and  wind- 
speed  (ms  ’)  from  the  Navy  Operational  Global  Atmospheric 
Prediction  System  (NOGAPS)  and  the  CCMP  at  51AE,  27°N  revealed 
that  the  winds  were  accurate  in  the  NOGAPS  forcing  (Fig.  2).  The 
frequent  changes  in  the  wind  direction  were  in  good  agreement 
with  the  CCMP.  This  suggests  that  the  model  responses  to  shamal 
winds  are  realistic  and  thereby  provided  confidence  in  our  model 
results.  The  shamal  winds  started  blowing  on  23  November  18Z 
(13ms'1),  reached  peak  strength  on  24  November  12Z 
(16ms  *)  and  dropped  to  2  m  s  1  on  27  November  00Z,  The 
winds  had  a  southerly  component  during  this  brief  period.  After 
this  short  break  the  wind-speed  picked  up  again  and  reached 
16  ms  1  on  28  November  00Z  followed  by  a  gradual  weakening. 


Nov  2004  Dec  2004 

Fig.  2.  Comparison  of  cross-calibrated,  multi-platform  (CCMP,  Atlas  et  al..  2009) 
(black)  and  NOGAPS  0.5  (grey)  wind-stress  (N  m  2)  and  wind-speed  (ms  *)  at 
51  E,  27N  (Fig.  la  for  location).  CCMP  wind  is  6  hourly  and  NOGAPS  is  3  hourly. 
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Fig.  1.  Surface  wind-speed  (shaded,  ms  ’)  and  vectors  from  the  cross-calibrated,  multi-platform  (CCMP.  Atlas  et  al..  2009)  before  and  during  a  winter  shamal  event 
(a)  20  November  and  (b)  12Z  24  November  2004. 
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By  4  December  OOZ  its  strength  decreased  to  6  m  s  \  It  should  be 
noted  that  the  CCMP  winds  are  six-hourly  and  the  NOCAPS  winds 
are  three-hourly. 

3.  Model  description 

The  HYCOM  (Bleck  et  al.,  2002)  is  used  in  order  to  simulate  the 
oceanic  responses  induced  by  the  shamal  winds.  Details  of  the 
model  configuration  and  forcing  fields  are  described  at  length  by 
Thoppii  and  Hogan  (2009)  and  therefore  only  briefly  presented 
here.  The  model  domain  extends  northward  from  22.3°N  to 
30.7°N  and  eastward  from  47.5°E  to  59.7°E.  The  model  has  16 
hybrid  vertical  layers  and  horizontal  resolution  of  ^1  km.  The 
baroclinic  (barotropic)  time-step  is  60s  (3s).  The  eastern 
boundary  is  closed  (59.7°E)  and  outfitted  with  a  50  grid  point 
buffer  zone  (~0.5°)  in  which  temperature,  salinity  and  pressure 
are  linearly  relaxed  toward  the  Generalized  Digital  Environmental 
Model  version  3.0  (GDEM3)  seasonally  varying  climatological 
values  (Newtonian  damping)  with  an  e-folding  time-scale  of  1  -76 
days.  The  monthly  river  inflow  into  the  basin  is  represented  as 
bogused  surface  precipitation.  The  model  treats  rivers  as  a  runoff 
addition  to  the  surface  precipitation  field,  which  acts  as  a 
freshwater  surface  flux  that  in  turn  decreases  the  salinity  of  the 
top  layer.  The  monthly  mean  river  discharge  values  are  specified 
at  the  river  mouths  (nearest  grid  cells).  Sea  surface  salinity  was 
restored  to  the  GDEM3  monthly  climatology  with  a  time-scale  of 
~30  days.  The  simulations  utilized  the  GISS  (NASA  Goddard 
Institute  for  Space  Studies  level  2)  vertical  mixing  scheme,  which 
is  constructed  using  the  Reynolds  stress  model  (Canuto  et  al., 
2001,  2002). 

The  model  integration  was  started  from  rest  on  January  1, 
2000  using  GDEM3  temperature  and  salinity  and  integrated  for 
six  years  for  the  period  2000-2005  using  1°,  three-hourly  NOGAPS 
forcing.  The  surface  forcing  includes  wind  speed,  vector  wind- 
stress,  air  temperature,  humidity,  precipitation,  surface  shortwave 
(daily  mean)  and  long-wave  heat  fluxes.  Surface  latent  and 
sensible  heat  fluxes,  along  with  evaporation  are  calculated  using 
bulk  formulae  (Kara  et  al.,  2002)  during  the  model  run  time  using 
model  SST.  The  model  reached  a  steady  seasonal  cycle  during  this 
period.  Furthermore,  the  model  integration  is  continued  for 
another  five  years  (2003-2007)  using  three-hourly,  1/2°  NOGAPS 
forcing  (which  is  available  only  after  2003)  in  which  the  winds 
were  corrected  for  the  land-sea  contamination  (Kara  et  al.,  2007) 
and  wind-speed  is  corrected  using  the  QuikScat  winds.  The  results 
presented  here  are  from  the  latter  run  using  1/2  NOGAPS  forcing 
unless  otherwise  indicated.  The  model  output  is  saved  as  daily 
snapshots  at  OOZ,  from  which  monthly  averages  are  prepared.  The 
sampling  once  per  day  is  not  likely  to  cause  major  bias  on 
the  mean  heat  flux  because  the  day-night  heat  flux  differences 
are  small  due  to  daily  mean  solar  radiation  being  used  (no  diurnal 
cycle).  The  results  presented  here  are  from  the  year  2004, 
particularly  during  the  November-December  2004  period. 

4.  Results 

4.1.  Seasonal  cycles  of  SST  and  heat  flux  during  2004 

Fig.  3  shows  the  seasonal  cycle  of  model  SST  and  heat  flux 
during  2004  and  comparisons  with  the  Moderate  Resolution 
Imaging  Spectroradiometer  (MODIS)  SST  and  heat  flux  based  on 
the  objectively  analyzed  air-sea  fluxes  (OAFIux,  Yu  and  Weller, 
2007).  The  agreement  between  the  model  and  MODIS  SSTs  are 
very  good  (correlation =0.99)  with  model  SSTs  being  cooler  by 
1  -2  "C  during  winter.  The  seasonal  cycle  in  the  SST  undergoes 


2004 

Fig.  3.  Seasonal  cycle  of  monthly  mean  SST  (  C)  from  model  (solid  dark  line)  and 
MODIS  (solid  grey  line)  and  net  surface  heat  flux  (W  m  2)  from  model  (dashed 
dark  line)  and  OAFIux  (dashed  grey  line)  (Yu  and  Weller,  2007)  in  the  Persian  Gulf 
region  (west  of  56  E)  during  2004.  The  basin  averaged  annual  mean  heat  flux  from 
the  model  (OAFIux)  is  -2.3  W  m  2  ( -  29  W  m  2), 


cooling  during  winter  (20  °C)  and  warming  in  summer  (32.5  °C) 
following  the  annual  cycle  of  the  net  surface  heat  flux  (NHF). 
Winter  cooling  peaks  in  February  which  is  two  months  after  the 
maximum  surface  heat  loss  of  200  Wm  2(145Wm  2 )  occurred 
in  the  model  (OAFIux)  in  December.  The  reduction  in  solar 
radiation  together  with  large  latent  heat  loss  (evaporation)  driven 
by  the  dry-air  emanating  from  the  northwest  results  in  surface 
heat  loss.  The  resulting  convective  vertical  mixing  decreases  the 
SST  below  ~20°C  in  the  Gulf  and  drives  a  deep  mixed  layer 
reaching  the  bottom  (not  shown).  During  summer  (May-August), 
the  upper-ocean  gains  heat  mainly  due  to  increased  solar 
radiation  concomitant  with  reduced  latent-heat  loss  owing  to 
weaker  winds.  The  upper-ocean  during  May-June  gains  heat  by 
108  W  m  2  which  is  higher  than  the  OAFIux  of  74  W  m  2.  The 
strong  near-surface  stratification  and  weak  wind-driven  mixing 
during  the  summer  produce  a  thin  mixed-layer  (~10m,  not 
shown)  and  warm  SSTs. 

It  is  clear  that  there  is  a  discrepancy  between  the  NHF  derived 
from  the  OAFIux  fields  and  that  obtained  from  the  model  (Fig.  3), 
especially  in  summer.  Despite  this  differences,  SSTs  during  this 
period  show  very  good  agreement.  The  basin-averaged  (west  of 
56  E)  annual  mean  NHF  from  the  model  is  -2.3  W  m“2.  This  value 
is  closer  to  the  basin-averaged  net-heat  loss  of  -7±4Wm  2 
obtained  from  the  estimates  of  advective  heat  fluxes  through  the 
Strait  of  Hormuz  by  Johns  et  al.  (2003)  and  the  adjusted  South¬ 
ampton  Oceanography  Centre  climatology  (SOC)  annual  mean  value 
of  4  W  m  2  (see  Johns  et  al.,  2003  for  details).  However,  the  basin- 
averaged  annual  mean  NHF  from  the  OAFIux  is  -29  W  m  2  which 
lies  outside  the  margin  of  error.  Caution  is  advised  when  interpret¬ 
ing  the  model-data  differences  because,  in  the  OAFIux,  there  are  only 
about  26,  1°  x  1°  grid  points  in  the  region  west  of  56°E 


42.  SST  and  heat-fluxes  during  a  shamal  event 

Time-series  of  the  basin  averaged  (west  of  56  E)  NHF,  latent 
and  sensible  heat  fluxes  from  the  OAFIux  during  a  shamal  event 
(November-December,  2004)  are  depicted  in  Fig.  4.  Pre-storm 
NHF  indicated  a  net  loss  of  90 Wm  2  heat  by  the  ocean  on 
November  22.  During  the  shamal  event,  the  Persian  Gulf 
experienced  strong  surface  heat  loss  of  474 Wm  2.  As 
expected,  latent  heat  dominated  the  total  heat  loss  and 
contribution  from  the  sensible  heat  flux  was  minimal  (Fig.  4a). 
In  order  to  understand  the  factors  contributing  to  the  total  heat 
flux,  we  examined  OAFIux  air-sea  temperature  differences  (Ts-Ta) 
and  specific  humidity  of  the  air  (qa,  g  kg  ’)  (Fig.  4b).  The 
northwesterly  winds  brought  cold  dry  air  into  the  Persian  Gulf 
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Fig.  4.  Time-series  of  basin  averaged  (a)  sensible  heat  flux  (SHF,  Wm  2),  latent  heat  flux  (LHF,  W  m  2),  net-heat  flux  (NHF,  Wm  2)  and  (b)  sea-air  temperature 
differences  {T,-!,*  C),  specific  humidity  of  air  ( qa ,  g  kg  ')  from  the  OAFlux  during  November- December  2004  period.  Model  NHF  (daily,  grey  bars)  and  NOGAPS  specific 
humidity  (3  hourly,  grey  line)  are  included  for  comparisons.  Note  that  the  OAFlux  Is  a  daily  mean  while  the  model  NHF  Is  a  daily  snapshot  at  00Z. 


region.  As  a  result,  the  specific  humidity  decreased  from 
14.7  g  kg  1  on  22  November  to  9.8  g  kg  1  on  25  November. 
This  sharp  drop  in  the  humidity  combined  with  strong  winds 
(16  ms  l)  enhanced  the  latent  heat  loss  (-433  Wm"2  on 
24  November).  A  comparison  of  specific  humidity  derived  from 
the  OAFlux  with  that  derived  from  the  1/2°  NOCAPS  revealed 
good  agreement  in  particular  a  rapid  drop  in  the  humidity 
associated  with  the  shamal.  The  air  temperature  also  dropped 
sharply  from  26.3  °C  on  22  November  to  20.7  °C  on  24  November 
resulting  in  6.3  °C  sea-air  temperature  difference.  This  led  to  an 
increase  in  the  sensible  heat  flux  of  104  W  m"2  during  the  storm 
(24  November).  Thus,  the  cold  dry  air  emanating  from  the 
northwest  enhanced  the  latent  heat  loss  which  led  to  a  net 
surface  heat  loss.  A  perturbation  experiment  with  a  constant 
Qa=  1 5  g  kg  1  reaffirmed  the  important  role  played  by  this  forcing 
parameter  on  the  total  heat  loss  (Section  4.6). 

The  surface  heat  flux  from  the  model  shown  in  Fig.  4a  (grey 
boxes)  qualitatively  agreed  with  the  OAFlux  (solid  line).  The 
model  NHF  was  larger  than  the  observations  especially  during 
the  shamal.  Part  of  the  difference  stemmed  from  the  fact  that  the 
model  NHF  is  a  snapshot  at  00Z  while  OAFlux  values  are  daily 
averages.  In  response  to  the  strong  wind  the  heat  loss  increased 
from  the  pre-storm  value  of  -27 Wm”2  (22  November)  to  a 
maximum  value  of  -671  W  m  2  during  25  November.  Both  the 
model  and  observations  exhibited  a  second  peak  of  heat  flux  on 
29  November  coinciding  with  a  period  of  increased  wind-speed 
(Fig.  2).  There  was  yet  another  period  of  enhanced  heat  loss 
which  occurred  around  17  December,  dominated  by  latent  heat 
loss,  however,  this  event  was  not  the  focus  of  this  study. 
Following  the  passage  of  the  shamal  winds  the  heat  loss  gradually 
decreased  to  a  near  pre-storm  condition.  The  mean  heat  flux  error 
between  the  model  and  the  observation  during  23-29  November 
was  -95  W  m  2. 

The  surface  heat  loss  resulting  from  the  shamal  winds  caused 
SST  cooling.  Figs.  5a  and  b  contrast  the  MOD1S  SST  (1  km)  a  few 
days  before  (17  November)  and  after  (29  November)  the  shamal 
event  (earliest  near  cloud-free  images  available)  and  their 
differences  (Fig.  5c)  quantified  the  SST  cooling.  The  SSTs  prior 
to  the  event  were  warmer  than  ~-28  C  over  most  of  the  region. 


The  shamal  event  caused  significant  SST  cooling  most  notably  in 
the  northwestern  and  shallower  shelf  regions.  The  SST  over  the 
shelf  was  noticeably  colder  by  up  to  6-8  C  due  to  wind-driven 
mixing  and  shallower  water  depth.  The  cool  dry  air  emanated 
from  the  land  enhanced  SST  cooling  below  17  X  in  the  northern¬ 
most  portion  of  the  Gulf  through  stronger  latent  heat  release. 
SST  cooling  of  2-4  °C  was  noted  in  the  northwestern  Gulf  where 
SST  dropped  below  25  °C.  Overall,  the  response  of  the  SST  to  the 
shamal  wind  was  stronger  in  the  northwestern  Gulf  and  was 
weaker  east  of  52°E  despite  strong  winds. 

The  model  reproduced  the  observed  SST  cooling  reasonably 
well.  Model  SST  difference  between  29  and  1 7  November  depicted 
in  Fig.  6a  were  in  accordance  with  the  MOD1S  SST  difference 
(Fig.  5c).  Maximum  SST  cooling  of  6-8  X  occurred  over  the  shelf 
with  magnitude  decreasing  offshore  (-4  to  -2  °C).  The  SST  near 
the  coast  was  colder  by  up  to  10  X,  which  could  not  be  confirmed 
in  the  MOD1S  due  to  cloud  cover.  Again,  the  northwestern  Gulf 
experienced  a  stronger  cooling  (2-4  X)  than  the  southeast 
(1-2  X).  This  differential  SST  cooling  between  the  northwestern 
and  southeastern  portions  of  the  Gulf  can  be  explained  in  terms  of 
contrasting  net-heat  flux  (Fig.  6b).  The  northwestern  region 
experienced  a  maximum  heat  loss  of  1400 Wm"2  while  the 
southeast  saw  only  a  400  W  m  2  heat  loss.  This  northwest  and 
southeast  contrast  in  heat  flux  can  be  attributed  to  the  gradient 
in  specific  humidity.  The  humidity  varied  from  8gkg  1  in  the 
northwest  to  13  g  kg  1  in  the  southeast  indicating  drier  air  and 
hence  larger  latent  heat  loss  in  the  northwest. 

To  quantify  the  atmospheric  forcing  that  led  to  the  SST  cooling, 
temperature  tendency  term  ( dT/dt  in  W  m  2)  and  total  heat  flux 
at  a  location  51°E,  27°N  are  presented  in  Fig.  7a.  As  expected,  SST 
changes  were  caused  in  large  part  by  the  surface  heat  flux.  The 
upper-ocean  experienced  a  net  heat  loss  in  excess  of  1000  W  m  2 
which  contributed  to  the  abrupt  SST  cooling.  The  net  heat  loss  and 
cool  SST  drove  a  deep  mixed-layer  of  25  m  (mixed-layer  increased 
from  36  to  61  m)  due  to  convective  vertical  mixing.  Note  that 
the  mixed-layer  reached  bottom  (30  November)  following  the 
passage  of  the  storm.  The  study  demonstrated  that  the  heat 
budget  balance  during  a  shamal  event  is  between  the  surface  heat 
loss  and  the  temperature  tendency  term.  The  basin  averaged 
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Fig.  5.  MODIS  (1  km)  derived  maps  of  SST  during  (a)  17  November  (b)  29 
November  2004  and  (c)  SST  differences  (A55T=SSTnoV29-SSTriovl7)  depicting  the 
SST  contrasts  before  and  after  the  shamal  event.  Black  areas  indicated  the  cloud 
where  data  were  not  available.  The  SST  over  the  shelf  was  noticeably  colder  by  up 
to  6-8  °C. 

(west  of  56JE)  contribution  of  net  heat  flux  to  the  SST  changes 
over  time  (time  integral)  is  shown  in  Fig.  7b  in  terms  of  the 
equivalent  cooling  in  °C.  This  emphasizes  that  the  water  column 
cooled  rapidly  from  the  air-sea  heat  loss  associated  with  the 
shamal.  The  cumulative  SST  changes  between  20  and  30 
November  were  -3.6  C.  The  SST  cooling  closely  followed  the 
surface  heat  loss,  suggesting  that  the  advective  and  diffusive 
processes  were  negligibly  small. 

4.3.  Vertical  sfrari/icarion 

The  shamal  winds  not  only  generated  SST  cooling  but  also 
caused  the  erosion  of  thermocline.  Vertical  profiles  of  tempera¬ 
ture  and  salinity  at  two  locations  in  the  Gulf  (54  E,  26.2  N  and 
5 I  E,  27  N),  for  three  days  (20,  27  November  and  5  December), 
shown  in  Fig.  8  illustrate  this.  Before  the  shamal  event,  the  water 


column  was  strongly  stratified  at  both  locations  with  a  well 
defined  mixed  layer  depth  of  30  m.  After  the  event,  the 
thermocline  at  the  northwestern  location  (51  °E.  27^N,  Fig.  8b) 
eroded  and  the  entire  water  column  mixed  down  to  the  bottom 
(60  m)  by  5  December.  Between  20  November  and  5  December, 
the  mixed-layer  temperature  decreased  by  3.5  °C.  Similarly  the 
salinity  stratification  also  underwent  changes  associated  with  the 
shamal  winds.  The  pre-storm  salinity  increased  with  depth. 
The  mixed-layer  salinity  increased  from  the  pre-shamal  event 
by  ~0.15%o  due  to  enhanced  evaporation  induced  by  the  shamal 
winds.  Contrarily,  the  shamal  had  a  lesser  impact  on  the 
stratification  at  the  southeast  location  (54°E,  26.2°N,  Fig.  8a) 
where  surface  heat  loss  was  weaker  (-200  Wm  2,  Fig.  6b).  The 
mixed-layer  deepened  and  cooled  by  2  °C.  The  stratification 
remained  stable  with  salinity  increasing  with  depth.  Thus,  the 
shamal  event  during  24  November  accelerated  winter  cooling  and 
convective  mixing  due  to  surface  heat  loss  driven  by  strong  winds 
and  cold  dry  air  emanated  from  the  northwest. 


4.4.  Circu/arion  and  Eddies 

One  of  the  important  features  of  circulation  in  the  Persian  Gulf 
is  the  boundary  intensified  northwestward  flowing  Iranian 
Coastal  Current  (ICC)  which  is  an  alongshore  flow  in  geostrophic 
equilibrium.  Between  April  and  July,  the  ICC  becomes  progres¬ 
sively  stronger  and  penetrates  far  into  the  northern  Gulf  due  to 
strengthening  stratification,  increasing  inflow  and  weakening 
winds.  The  ICC  evolves  into  a  series  of  mostly  cyclonic  eddies 
between  July  and  August.  A  total  of  three  to  five  eddies  exist  and 
their  presence  dominate  the  circulation  during  the  August- 
September  period  (Thoppil  and  Hogan,  submitted  for  publication). 
Because  of  the  narrow,  shallow  and  irregular  nature  of  bathy¬ 
metry,  these  eddies  remain  stationary  or  trapped  by  the  bottom 
topography  until  they  dissipate  locally.  They  have  clear  signatures 
in  the  surface  currents  until  November,  when  winter  cooling 
causes  the  erosion  of  stratification  and  the  collapse  of  eddies.  The 
year  2004  was  not  an  exception;  the  remnants  of  eddies  could  be 
seen  prior  to  the  shamal.  Fig.  9  shows  the  MODIS  (1  km)  SST, 
chlorophyll  concentration  (mgm  3)  and  the  model  mixed-layer 
currents  during  10  November  2004.  The  model  currents  indicated 
two  cyclonic  eddies  (52.7  E,  26.8  N;  54.2°E,  26.1°N)  and  an 
anticyclonic  eddy  (53.5°E,  26.2  3N)  sandwiched  between  them 
with  their  diameters  ranging  from  ~50to  100  km.  The  model  also 
shows  a  string  of  small  cyclonic  and  anticyclonic  eddies  traveling 
from  the  Strait  of  Hormuz  into  the  Persian  Gulf.  The  signature  of 
two  cyclonic  eddies  can  be  delineated  from  the  images  of  MODIS 
SST  and  chlorophyll  observations  while  smaller  eddies  were 
difficult  to  detect. 

The  strong  winds  associated  with  the  shamal  event  had  a 
prominent  signature  in  the  current  field.  The  mixed-layer  currents 
prior  to  the  shamal  event  (20  November)  resembled  those  during 
10  November  (Fig.  10)  in  that  the  remnants  of  cyclonic  eddies 
(from  September  2004)  were  the  dominant  circulation  features. 
There  was  no  indication  of  a  northwestward  flowing  ICC  at  this 
time.  As  expected  the  strong  shamal  winds  intensified  the 
circulation  in  the  entire  Gulf,  while  the  overall  responses  were 
different  in  the  northern  and  southern  (alongshore  direction) 
portions  of  the  Gulf.  The  shamal  established  an  ICC  with  a  distinct 
maximum  offshore  speed  (25  cm  s  1 ),  from  the  Strait  of  Hormuz 
to  about  52°E,  where  it  veered  offshore.  While  ICC  is  a  persistent 
feature  between  April  and  July  (Thoppil  and  Hogan,  submitted  for 
publication)  its  reappearance  during  the  shamal  was  indeed 
interesting.  We  also  note  the  strengthening  of  the  southeastward 
current  along  the  southern  portion  from  Saudi  Arabia  to  the 
United  Arab  Emirates,  obviously  driven  by  shamal  winds.  A  partial 
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Fig.  6.  Model  (a)  SST  differences  (°C)  between  17  November  and  29  November  2004  (ASST =SSTnov2Q-SSTmyVn)  and  (b)  net  heat  flux  (W  m  2)  dunng  24  November  2004 
from  the  control  run  (cr).  (c-d)  Same  as  above  but  from  the  perturbation  experiment  (pe)  in  which  specific  humidity  was  specified  to  have  a  constant  value  of 
qa  =  \5  g  kg  \  The  contours  are  for  (c)  SST  cooling  differences  between  the  two  experiments  (ASST^-ASST^)  with  an  interval  of  1  X  (grey  contour  being  the  1  X)  (d)  NHF 
differences  between  the  two  experiments  (NHF^-NHF^)  during  24  November  2004  with  an  interval  of  100  Wm  2  (grey  contour  being  the  500  W  m  2).  That  is,  in  the 
constant  humidity  run,  SST  cooling  was  reduced  by  1  X  and  net  heat  loss  was  reduced  by  500  Wm  2  from  the  control  run. 
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Fig.  7.  (a)  Rate  of  change  of  SST  (?T/?f.  Wm  2).  net-heat  flux  (NHF,  W  m  2)and  mixed-layer  depth  (MLD,  m)  at  SI  E,  27  N  (Fig.  la  for  location)  and  (b)  the  basin  averaged 
(west  of  56  E)  time  integral  (cumulative)  of  SST  change  and  net-heat  flux  expressed  in  X.  The  SST  cooling  closely  followed  the  air-sea  heat  loss. 


relaxation  of  winds  on  26  November  (Fig.  2)  caused  the  currents 
to  return  to  a  near  pre-storm  condition  on  27  November.  The 
winds  continued  forcing  a  persistent  ICC  until  6  December.  On 
5  December  the  ICC  and  a  returning  flow  to  the  south  constituted 
a  closed  cyclonic  circulation  cell  between  52°E  and  54  E  and  this 
cell  broke  down  into  smaller  eddies  during  7  December  (Fig.  10). 
Summarizing  the  shamal  response,  currents  along  the  southern 
portions  (alongshore  direction)  of  the  Gulf  were  directly  wind- 
driven  while  the  currents  along  the  northern  portions  were  the 
indirect  result  of  the  shamal  winds  through  a  lagged  response  to 
the  sea  surface  height  (SSH)  difference. 


4.5.  Sea  surface  height  and  transport 

We  now  examine  the  basin  wide  barotropic  response  to  the 
strong  winds.  The  simulated  SSH  before  and  after  the  shamal 
event  is  depicted  in  Figs.  1  la  and  b.  As  a  consequence  of  the  wind 
the  SSH  minimum  occurred  in  the  northernmost  part  of  the  Gulf 
where  SSH  dropped  below  20-30  cm.  On  the  other  hand,  SSH 
increased  by  10-20  cm  in  the  eastern  half  of  the  basin.  The 
strongest  SSH  slope  appeared  in  the  northernmost  and  southern 
shallow  regions.  In  the  latter  case,  the  water  piled  up  along  the 
coast  in  the  southern  shallow  regions  by  a  strong  wind  formed  a 
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Fig.  8.  Vertical  profiles  of  temperature  (solid  lines)  and  salinity  (dashed  lines)  at  two  locations;  54  E,  26.2  N  (left)  and  51  E,  27  N  (right)  during  20  November, 
27  November  and  5  December  2004.  See  Fig.  la  for  locations. 
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Fig.  9.  MODIS  (1  km)  derived  SST  (X).  chlorophyll  concentration  (mgm  3)  and 
model  surface  currents  overlain  on  the  model  SST  during  10  November  2004. 
Ocean  color  provides  a  better  depiction  of  eddies  and  their  location  compared  well 
with  that  in  the  model.  It  suffices  to  mention  that  green  represents  lower 
concentrations  and  yellow  and  red  higher  concentrations,  since  we  are  not 
concerned  with  the  magnitude  of  the  chlorophyll  concentrations. 


high  sea  surface  of  20-30  cm.  The  SSH  pattern  simulated  by  the 
model  qualitatively  agreed  with  a  storm  surge  model  responding 
to  an  idealized  shamal  wind  structure  (El-Sabh  and  Murty,  1989). 
They  noted  positive  surges  up  to  2-3  m  in  the  shallow  southern 
regions  and  Bay  of  Bahrain  and  negative  surges  of  1  m  in  the 
northernmost  gulf  (their  Fig.  11),  although  it  should  be  noted  that 
their  model  included  tidal  forcing. 

Sea  level  observation  from  a  coastal  station  near  Bahrain 
(50.6°E,  26.2  N)  was  used  to  examine  the  coastal  response  to  the 


shamal  wind.  The  observed  sea  level  anomaly  (SLA,  after 
subtracting  the  annual  mean  of  2004)  is  plotted  (grey  bars)  along 
with  model  SSH  at  51  °E,  27  N  (solid  line)  in  Fig.  11c.  The  SLA  was 
positive  prior  to  the  storm  (~10cm).  The  shamal  wind  caused  a 
drop  of  —16  cm  in  SLA  on  25  November  followed  by  a  rise  of 
18  cm  in  SLA  on  27  November.  The  model  SSH  at  this  location 
failed  to  simulate  the  observed  sea  level  pattern.  Instead,  an 
examination  of  SSH  farther  north  (51  °E,  27°N)  revealed  a 
qualitative  representation  of  the  observed  sea  level  variability. 
The  fluctuations  of  the  sea  level  were  in  phase  with  the  model 
SSH,  an  indication  that  the  sea  levels  were  being  influenced  by 
changes  in  wind  forcing  (Fig.  2).  It  should  be  noted  that  tidal 
components  are  also  important  in  determining  the  total  sea  level 
in  the  Persian  Gulf,  which  has  not  included  in  the  model. 

The  sea  level  difference  between  the  northwestern  (west  of 
52 °E)  and  southeastern  (52E-56  E)  portions  of  the  Gulf  induced 
by  shamal  winds  can  be,  in  part,  responsible  for  the  variability  in 
the  ICC.  The  area  averaged  SSH  difference  (ASSH=SSHwrsf-SSHeasf) 
and  zonal  velocity  (u)  at  53.1°E,  26.7°N  within  the  core  of  ICC 
depicted  in  Fig.  12a  showed  significant  correlation  of  0.72  with  u 
leading  AS5H  by  1  day.  It  is  also  reasonable  to  speculate  that  the 
flow  variability  can  arise  from  the  variability  in  the  inflow- 
outflow  exchange  through  the  Strait  of  Hormuz.  The  transport 
through  the  Strait  at  56.2  E  for  the  northern  (north  of  26.4nN)  and 
southern  portions  (south  of  26.4°N)  clearly  showed  variability 
associated  with  the  shamal  (Fig.  12b).  The  outflow  increased  to 
0.25  Sv  (Sv=  lO^m’s  *)  on  24  November  followed  by  an 
enhanced  inflow  of  0.25  Sv  during  25  November.  Similarly,  an 
accelerated  outflow  (0.17  Sv)  followed  by  an  inflow  (0.26  Sv) 
occurred  during  28-29  November.  These  periods  of  enhanced 
inflow  produced  a  stronger  ICC  along  the  coast  of  Iran  (Fig.  10) 
implying  that  they  are  connected.  We  speculate  that  the  ICC  is  a 
lagged  response  to  the  SSH  difference,  however,  a  detailed 
analysis  is  beyond  the  scope  of  this  paper. 


4.6.  Sensitivity  to  humidity 

In  order  to  identify  and  isolate  the  influence  of  humidity  on  the 
latent  heat  loss  and  hence  on  the  net  heat  loss  during  a  shamal 
(Figs.  4  and  6b),  we  carried  out  a  perturbation  model  experiment 
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Fig.  10.  Model  mixed-layer  currents  superimposed  on  the  SST  during  20  November,  25  November,  27  November.  29  November,  3  December  and  7  December  2004.  For 
clarity,  velocity  vectors  whose  magnitudes  larger  than  4  cm  s  1  are  plotted.  Every  10th  vector  is  plotted.  Iranian  Coastal  Current  (ICC)  is  marked. 
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Fig.  11.  Sea-surface  height  (SSH,  cm)  during  (a)  20  November  and  (b)  25  November  2004.  The  contour  interval  is  5  cm  and  the  thick-line  denotes  zero  contour,  (c) 
Observations  of  daily  mean  sea  level  anomaly  (anomaly  is  obtained  by  subtracting  the  annual  mean  of  2004)  from  a  coastal  station  near  Bahrain  (50.6  F,  26.2  N;  grey  bars) 
and  model  SSH  at  5 I  E.  27  N  (solid  line;  see  Fig.  1 1b  for  locations). 


by  maintaining  a  constant  value  of  qa  =  15gkg  1  specific 
humidity  of  air  from  1  to  30  November  2004.  The  choice  of  this 
value  broadly  represented  the  pre-storm  conditions  in  both 
OAFlux  and  the  NOGAPS  fields  (Fig.  4b).  The  SST  differences 
(29-17  November)  and  net  heat  flux  (24  November)  from  this 
experiment  are  depicted  in  Figs.  6c  and  d.  The  difference  in  net 
heat  flux  in  the  constant  humidity  run  was  substantial  and 
accounted  for  a  change  of  400-500  W  m~2  from  the  control  run. 


The  northwestern  Gulf  experienced  a  maximum  net  heat  loss  of 
800  Wm"2  in  contrast  to  1400  W  m  2  in  the  control  run.  These 
differences  were  primarily  due  to  a  reduction  in  the  latent  heat 
flux  caused  by  higher  humidity.  When  averaged  over  the  basin 
west  of  56  E,  the  NHF  due  to  humidity  alone  accounted  for  38%  of 
the  total  heat  loss.  The  reduction  of  NHF  had  a  noticeable  impact 
on  the  SST  cooling  especially  along  the  coastal  regions  (Fig.  6c), 
where  SST  cooling  was  weaker  than  the  control  run  by  2-4  °C. 
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Fig.  12.  (a)  Sea  surface  height  difference  between  the  northwestern  (west  of  52"E) 
and  southeastern  (52-56°E)  portions  of  the  Persian  Gulf  (dotted  line)  and  zonal 
velocity  (u,  cm  s  ')  within  the  core  of  Iranian  Coastal  Current  (ICC)  at  S3.1  E. 
26.7  N.  Note  that  SSH  difference  is  plotted  with  a  1-day  lag.  (b)  Transport  (Sv) 
across  the  Strait  of  Hormuz  at  56.2  E  for  the  northern  (north  of  26.4  N;  dark  bars) 
and  southern  portions  (south  of  26.4  N:  grey  bars)  of  the  Strait. 


Elsewhere  the  SST  cooling  was  0-1  X  different  from  the  control 
run.  Thus,  the  sensitivity  experiment  reiterated  the  role  of 
humidity  in  determining  the  winter  cooling  and  heat-budget  in 
the  northern  Gulf  during  a  shamal  event. 


5.  Conclusions 

Results  from  a  high-resolution  HYbrid  Coordinate  Ocean 
Model  (HYCOM),  forced  by  1/2°  NOGAPS  atmospheric  data,  were 
used  in  order  to  study  the  dynamical  response  of  the  Persian  Gulf 
to  a  winter  time  shamal  forcing.  The  period  of  20  November-5 
December  2004  was  selected  because  it  included  a  well  defined 
shamal  event  that  lasted  about  4-5  days  (24-29  November).  In 
addition  to  the  strong  winds  (16  ms-1),  the  winter  shamal 
brought  cold  dry  air  that  lead  to  enhanced  surface  heat  loss  and 
evaporation.  The  Persian  Gulf  experienced  strong  cooling  with 
surface  heat  loss  in  excess  of  lOOOWm  2.  This  resulted  in 
significant  SST  cooling  most  notably  in  the  northern  and 
shallower  shelf  regions.  The  SST  over  the  shelf  was  noticeably 
colder  by  up  to  10  X  due  to  surface  heat  loss  and  wind-driven 
mixing.  The  heat-budget  balance  during  the  shamal  event  was 
between  the  surface  heat  flux  and  the  temperature  tendency  term 
and  contributions  from  the  advective/diffusive  term  appeared 
negligible.  A  sensitivity  experiment  with  a  constant  specific 
humidity  of  ga=15gkg  1  confirmed  that  about  38%  of  net  heat 


loss  was  due  to  the  air-sea  humidity  differences.  The  shamal 
event  also  significantly  contributed  to  the  erosion  of  vertical 
stratification.  The  pre-shamal  vertical  temperature  at  51  "E,  27nN 
indicated  a  strongly  stratified  water  column  with  a  mixed  layer 
depth  of  30  m.  After  the  shamal  thermocline  eroded  and  the 
entire  water  column  mixed  down  to  its  bottom  (60  m). 

The  remnants  of  mesoscale  to  sub-mesoscale  eddies  from  late 
summer  (August-September  2004)  dominated  pre-shamal  circu¬ 
lation  in  the  Gulf.  Perhaps,  one  of  the  most  intriguing  features  of 
the  circulation  during  the  shamal  was  the  genesis  of  a  25  cm  s  1 
northwestward  flowing  Iranian  Coastal  Current  (ICC)  from  the 
Strait  of  Hormuz  to  north  of  Qatar.  It  veered  offshore  at  about 
52X  with  a  branch  continuing  northwestward  and  another 
branch  joining  the  return  southeastward  flow  to  the  right. 
The  shamal  winds  also  intensified  the  outflowing  branch  of  the 
cyclonic  circulation;  part  of  it  terminated  before  reaching  the 
Strait.  The  ICC  persisted  until  6  December  and  then  collapsed  into 
smaller  eddies  due  to  weakening  winds  and  eroding  stratification. 

The  barotropic  response  of  the  Gulf  to  shamal  event  reasonably 
agreed  with  the  observations  of  sea  level.  The  strongest  positive 
sea  level  slope  was  generated  in  the  southern  shallow  regions 
associated  with  the  wind  setup  against  the  United  Arab  Emirates 
coast  (piling  up  of  water  near  the  coast)  while  the  strongest 
negative  sea  level  slope  was  formed  in  the  northernmost  part  of 
the  Gulf.  The  area  averaged  sea  surface  height  (SSH)  difference 
and  zonal  velocity  (u)  within  the  core  of  ICC  showed  a  significant 
correlation  of  0.72  with  u  leading  SSH  by  1-day  suggesting  that 
part  of  the  ICC  variability  was  driven  by  the  SSH  fluctuations. 

We  have  for  the  first  time  presented  the  dynamic  response  of 
the  Persian  Gulf  to  a  wintertime  shamal  event  using  model  results 
and  satellite  observations  with  a  special  emphasis  on  the  surface 
circulation  and  heat-budget.  The  shamal  winds  in  summer  will 
induce  a  different  response  because  summertime  stratification  is 
stronger  and  has  a  much  stronger  inflow  of  low-salinity  water 
along  the  coast  of  Iran.  The  extent  to  which  the  short-lived  shamal 
winds  impact  the  long-term  circulation  in  the  Gulf  is  relatively 
unknown.  The  sharp  differences  in  the  stratification  between 
summer  and  winter  make  the  response  time  differ.  Although  the 
prevailing  winds  are  generally  northwesterly  throughout  the  year, 
short  episodes  of  very  strong  shamal  winds  can  change  the 
surface  circulation  abruptly  during  short  spells  of  a  few  days  to 
weeks.  Frequent  episodes  of  shamal  winds  may  lead  to  an 
increase  in  the  intensity  of  the  cyclonic  circulation  by  altering 
the  exchange  circulation  through  the  Strait  of  Hormuz.  Exclusion 
of  these  events  in  the  modeling  studies  are  likely  to  under¬ 
estimate  the  surface  circulation.  More  systematic  studies  using 
remotely  sensed  and  in  situ  data  in  conjunction  with  numerical 
modeling  would  aid  in  our  understanding  of  the  complex  air-sea 
interaction  and  circulation  that  are  taking  place  during  shamal 
events.  A  prior  knowledge  of  circulation  in  response  to  the  shamal 
wind  is  not  only  useful  for  navigation,  but  is  essential  to 
understanding  of  the  movement  of  oil  slicks. 
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